Dynamics and transport of the spin liquid: 
application to k-(ET)2Cu2(CN)3 
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We describe neutron scattering, NMR relaxation, and thermal transport properties of Z2 spin 
liquids in two dimensions. Comparison to recent experiments on the spin S = 1/2 triangular lattice 
antiferromagnet in k-(ET)2Cu2(CN)3 shows that this compound may realize a Z2 spin liquid. We 
argue that the topological 'vison' excitations dominate thermal transport, and that recent thermal 
conductivity experiments by M. Yamashita et al. have observed the vison gap. 



Much attention P, d d, H d, d, 13 has recently focused 
on the organic compound k-(ET)2Cu2(CN)3 because it 
may be the first experimental realization of a resonating 
valence bond spin liquid [1, Q ■ This compound belongs 
to a class [13, Ell of organic Mott insulators which can 
be described hy S — 1/2 spins residing on the vertices 
of a triangular lattice. Experiments have not detected 
any magnetic order or a structural distortion leading to 
a doubling of the unit cell in k-(ET)2Cu2(CN)3, and so 
there is justifiable optimism that the elusive spin liquid 
state may finally have been found. 

The debate then turns to the identification of the pre- 
cise spin liquid state, among the many possible candi- 
dates. Measurements of the electronic specific heat, Cp, 
by S. Yamashita et al. Q were interpreted to yield a non- 
zero low temperature (T) value of 7 = ImiT^o Cp/T. 
Such a non-zero 7 is characteristic of a Fermi surface, 
and hence a spin liquid state with a Fermi surface of neu- 
tral, S — 1/2, fcrmionic spinous was postulated 0, d, d|- 
However, it should be noted that the measurement of 7 
involves a potentially dangerous subtraction of a diver- 
gent nuclear specific heat [3|. 

Very recently, M. Yamashita et al. have measured [ij 
the thermal conductivity, k, to below T « 0.1 K. This 
has the advantage of focusing on the mobile excitations, 
and not being contaminated by a nuclear contribution. 
A spinon Fermi surface should yield a non-zero low T 
limit for k/T, but this quantity was clearly observed to 
vanish. Instead, the measured k was fit reasonably well 
by the activated behavior k ~ exp(— A^/T), with a 'gap' 
Ak « 0.46 K. Furthermore, k was found to be insensitive 
to an applied field for iJ < 4 T, suggesting that the gap 
Ak is associated with a spinless excitation. These ob- 
servations appear to be incompatible with spinon Fermi 
surface states at these low T, and we shall present an 
alternative theory here. 

Also of interest are the measurements Q of the NMR 
relaxation rate, l/Ti. The power-law behavior 1/Ti ^ 
T°, with the exponent a ~ 1.5, was observed for 0.02 < 
T < 0.3 K. This requires the presence of spinful excita- 
tions with a gaplcss spectrum at the fields of the NMR 
experiment, although at zero field there may well be a 
small spin gap. 



In this paper, we will compare these observations 
with the Z2 smn liquid state originally proposed in 
Refs. llSl, llJ- The low energy excitations of this 
state are described by a Z2 gauge theory, and the spin- 
ful excitations are constructed from S — 1/2 quanta 
(the spinous) which carry a Z2 electric charge. Crucial 
to our purposes here are vortex-like spinless excitations 
15 1 which carry Z2 magnetic flux, later dubbed 'visons' 
Iq . A number of solvable models of Z2 spin liquids, 
with spinon and vison excitations, have been constructed 
S [13, H, [US mi [ii]. We propose here that it is 
the visons which dominate the thermal transport in k- 
(ET)2Cu2(CN)3, and the gap A„ is therefore identified 
with a vison energy gap, A„. If our interpretation is cor- 
rect, the vison has been observed by M. Yamashita et al. 

Our proposal requires that the density of states of 
low energy vison excitations is much larger than that 
of all other excitations. A model appropriate to k- 
(ET)2Cu2(CN)3 is the triangular lattice S" = 1/2 antifer- 
romagnet with nearest neighbor two-spin exchange (J2) 
and plaquette four-spin ( J4) exchange which was studied 
by Liming et al. [23|. They found antiferromagnetic or- 
der at J4 = (as in earlier work (31), and a quantum 
phase transition to a spin liquid state with a spin gap 
around J4/J2 ~ 0.1. Notably, they found a very large 
density of low- lying spin singlet excitations near the tran- 
sition. We propose here that k-(ET)2Cu2(CN)3 is near 
this quantum phase transition, and identify these singlets 
with visons which have a small gap and bandwidth, both 
much smaller than the spin exchange J2 ~ 250K. We will 
argue below that at T ^ J2 , and comparable to the vison 
bandwidth, visons will dominate the thermal transport. 

Further support for the proximity of a magnetic order- 
ing quantum critical point comes from 11] the closely re- 
lated series of compounds X[Pd(dmit)2]2- By varying the 
anisotropy of the triangular lattice by varying X, we ob- 
tain compounds with decreasing magnetic ordering criti- 
cal temperatures, until we eventually reach a compound 
with a spin gap and valence bond solid (VBS) order ^^23\ . 
In between is the compound [27] with X=EtMe3P which 
has been proposed to be at the quantum critical point 
111 ], and has properties similar to k-(ET)2Cu2(CN)3. Fi- 
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nally, series expansion studies [26| also place the triangu- 
lar lattice antiferromagnet near a quantum critical point 
between magnetically ordered and VBS states. 

A description of the NMR experiments requires a the- 
ory for the spinon excitations of the Zi spin liquid. The 
many models of spin liquids 12 , [ij, [3, 15 , [3, 17 



critical form 



18, l2ll . |22| | have cases with either fermionic or 

bosonic spinous. While we do not find a satisfactory ex- 
planation for the NMR with fermionic spinous, we show 
that a model 12, [l^ id] of bosonic spinous in a spin 



liquid close to the quantum phase transition to the anti- 
ferromagnetically ordered state (as found in the model of 
Liming et al. (23[) does naturally explain the T depen- 
dence of 1/Ti. We shall show below that the quantum 
critical region for this transition leads to 1/Ti with 
the exponent [H, [2§| fj = 1.37, reasonably close to the 
measured value a = 1.5. It is important to note that the 
vison gap, A^, remains non-zero across this magnetic or- 
dering critical point [35|. Consequently, our interpreta- 
tion of the experiments remains valid even if the system 
acquires a small antiferromagnetic moment, as may be 
the case in the presence of the applied magnetic field 
present in the NMR measurements. 

The remainder of the paper presents a number of com- 
putations of the physical properties of Z2 spin liquids, 
and uses them to elaborate on the experimental inter- 
pretation sketched above. 

We begin with a theory [s^l of the spinon excitations 
near the quantum critical point between the magnetically 
ordered state and the Z2 spin liquid. Here the low energy 
spinous are 5 = 1/2 complex bosons z^, with a =tii ^ 
spin index, and the low energy imaginary time action is 



(1) 



where (r, r) are spacetime co-ordinates, 17 is a coupling 
which tunes the transition to the spin liquid present for 
some 5 > 5cj and c is a spin- wave velocity. We impose 
the local constraint J2a I-^qP = 1 ^i^^ of * quartic 
self- interaction between the spinous. This theory has an 
emergent 0(4) global symmetry 2^, sH (which becomes 
manifest when Za is written in terms of its real and imag- 
inary components). This symmetry is an enlargement of 
the SU(2) spin rotation symmetry, and we will neglect the 
irrelevant terms which reduce the symmetry to SU(2). 

(i) Dynamic spin susceptibility. The dynamic spin 
correlations of S near the quantum critical point can be 
computed by the 1/A^ expansion on the 0{N) model, 
which has been described elsewhere 32| . With an eye to- 
wards possible future neutron scattering measurements, 
we first describe the dynamic spin susceptibility, xikju) 
as a function of momentum k and real frequency co. Here 
the momentum k is measured as a deviation from the 
ordering wavevector, Q, of the antiferromagnetically or- 
dered state. At g = gc and T = 0, this has the quantum- 
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(2) 



where the exponent fj is related to the scaling dimension 
of the composite spin operator ZaCr^^cTjisZis {a are the 
Pauli matrices), and is known with high precision from 
field-theoretic studies [2^ {fj = 1.374(12)) and Monte 
Carlo simulations 2^ {fj = 1.373(2)). The overall am- 
plitude A is non-universal, but the same A will appear 
in a number of results below. Integrating Eq. ^ over 
all k, we obtain the local susceptibility Xi(cj), which is 
also often measured in scattering experiments, again at 
g = gc and T = 



ImxL(w) 



Asgn{Lu) sin(7r77/2) 



4c2 



7Tfj/2 



(3) 



Let us now move into the spin liquid state, with g > g^ 
where the spinous have an energy gap A^. The critical 
results in Eqs. ([2]) and ^ will apply for |a;| » Az, but 
for |ti;| ~ 2Az, we will have spectra characteristic of the 
creation of a pair of spinous (we set h = 1, although 
it appears explicitly in a few expressions below). Com- 
puting the pair creation amplitude of non-interacting 
spinous, we obtain a step-discontinuity threshold at lu — 
v/c2pT4Af (at T = 0). However, the spinons do have 
a repulsive interaction with each other, and this reduces 
the phase space for spinon creation at low momentum, 
as described in the supplementary material; the actual 
threshold behavior is: 



Imx(fc,cj) 



ACsgii{Lu) 



(1^1 -V^ 



4A5 



In^ 



UJ 



(4) 



16A2 



where C is a universal constant; to leading order in the 
1/A^ expansion, C — N'^ /16. We can also integrate the k- 
dependent generalization of Eq. ^ to obtain a threshold 
behavior for the local susceptibility at 2Az: lmxL{^) ~ 
sgn(cc;)(|a;|-2A,)/ln'(|cc;|-2A,). 

(a) NMR relaxation. Turning to the NMR relaxation 
rate, we have to consider T > 0, and compute 



r = lim — — lmxL{i^)- 



(5) 



This is far more subtle than the computations at T = 0, 
because we have to compute the damping of the quantum 
critical excitations at T > and extend to the regime 
u! ^ T. From general scaling arguments 34], we have 



F = :^{kBTfHAz/{kBT)), 
c^ 



(6) 



where $ is a universal function. The computation of 
$ for undamped spinons at iV = 00 is straightforward. 
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and unlike the case for confining antiferromagnets [32| . 
yields a reasonable non-zero answer: 
V4 + ey)]~^. However, the corrections are singular, 
because F has a singular dependence upon the spinon 
lifetime. A self-consistent treatment of the spinon damp- 
ing is described in the supplementary material, and leads 
to the quantum-critical result (A^ = 0): 



$(0) 



(V5-1) 
IGtt 



1 + 0.931- 



IniV 



N 



(7) 



(in) Thermal conductivity. We now turn to the ther- 
mal transport co-efficient measured in the recent reveal- 
ing experiments of Ref. '4. We consider the contribution 
of the spinous and visons in turn below, presenting fur- 
ther arguments on why the vison contribution can dom- 
inate in the experiments. 

(Hi. a) Spinous. For agreement with the NMR mea- 
surements of 1/Ti we need the spinous to be in the 
quantum critical regime, as described above. Therefore, 
we limit our considerations here to the quantum critical 
thermal conductivity of the spinous, k^, with = 0. 
This can be obtained from the recent general theory of 
quantum critical transport [33] which yields 



(8) 

where s is the entropy density of the spinous, and l/r^ 
is the spinon momentum relaxation rate, with the T de- 
pendence 



imp 



y2/i/- 



(9) 



Here i' is the critical exponent of the 0(4) model [34| . 
ly = 0.749(2), and so r™P ~ T-°-^^. The two dimen- 
sional entropy density can be obtained from the results 
of Ref. in 



3A^C(3)fc|T' 



0.3344 

N 



(10) 



where ^ is the Reimann zeta function. We estimate the 
co-efficient in Eq. ^ in the supplementary material using 
a soft-spin theory with the spinous moving in a random 
potential, y(r)|zQp, due to impurities of density njmp 
each exerting a Yukawa potential Vq = Vz/ [q^ + M^); this 
leads to 



2/iy-3 



(11) 



Here a is the spacing between the layers, and is the 
spinon bandwidth in temperature units and is propor- 
tional to the spinon velocity c. 

(iii.b) Visons. The visons are thermally excited 
across an energy gap, A^, and so can be considered to be 
a dilute Boltzmann gas of particles of mass m„ . We as- 
sume there are species of visons. The visons see the 



8 10 

I IT 



In(ft) 



FIG. 1: Fit of the T dependence of the vison thermal conduc- 
tivity in Eq. p2p to the thermal conductivity measurements 
by Yamashita et al. Tv, Av and the overall prefactor were 
the fit parameters. 



background filling of spins as a magnetic flux through 
the plaquette on the dual lattice, and hence the dynam- 
ics of visons can be well described by a fully-frustrated 
quantum Ising model on the honeycomb lattice. De- 
tailed calculations show that there are four minima of 
the vison band with an emergent 0(4) flavor symme- 
try at low energy [13], therefore Ny = 4. As with 
the spinous, the visons are assumed to scatter off im- 
purities of density rtimp with, say, a Yukawa potential 
Vq = Vv/{q^ + Ai^). We use the fact that at low T, 
and for a large vison mass m«, the visons are slowly 
moving. So each impurity scattering event can be de- 
scribed by a T- matrix = [rriy ln(l/fc)/7r]~^ characteristic 
of low momentum scattering in two dimensions. Applica- 
tion of Fermi's golden rule then yields a vison scattering 
rate 1/t™p = ir'^nin^p/ {my\n^{l/k)). This formula be- 
comes applicable when ln(l/A:) x Vy / {h^ fjL^ /2my) ^ 1 
i.e. the impurity potential becomes nonperturbative. 
We can now insert this scattering rate into a standard 
Boltzmann equation computation of the thermal con- 
ductivity Ky — 2kgTnyT™^^ /niy, whcrc Uy is the ther- 
mally excited vison density and the typical momentum 
k ^ {m.ykBTy^'^ , to obtain 



NyUiyklT^ \n{Ty/T)e 



-AJ(kBT) 



(12) 



Here Ty is some ultraviolet cutoff temperature which can 
be taken as the vison bandwidth. Note that for a large 
density of states of vison excitations, i.e. a large m„, the 
prefactor of the exponential can be large. Similar calcu- 
lations will not lead to a logarithmic divergence for the 
critical spinon z due to the positive anomalous dimension 
of ]z]2, and therefore the impurity scattering of spinous 
is perturbative for Vz/{cfj,h)'^ < 1. 

Using Eq. p^ . we fit the thermal conductivity mea- 
sured by M. Yamashita et al. in Ref. 3 by tuning param- 
eters Ty and A^. The best fit values are Ty = 8.15K, 
and Ay = — 0.238iir, as shown in Fig. [TJ For consis- 
tency check, we calculate the ratio between the thermal 
conductivities contributed by spinous and visons using 
Eq. (jlip and Eq. ([T^ and assuming moderate spinon 
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FIG. 2: Ratio of the thermal conductivity of spinons to visons 
in Eq. ([13} 



a small energy gap and bandwidth) dominate the thermal 
transport [J]. 

We are very grateful to Minora Yamashita for valuable 
discussions of the results of Ref. 0, and to the authors of 
Ref. for permission to use their data in Fig. [TJ We 
thank K. Kanoda, S. Kivelson, and T. Senthil for useful 
discussions. This research was supported by the NSF 
under grant DMR-0757145. 



impurity strength Vz/{cfih)'^ 



2/1/-3 



1 



1 



„A„/feBT 



(lnT„/r); 



We plot this ratio in Fig. [H with ^ J2 ^ 250 K 
and other parameters as above, for the experimentally 
relevant temperature between 0.1 K and 0.6 K; we find 
consistency because n is dominated by the vison contri- 
bution. The vison dispersion is quadratic above the vison 
gap, and this leads to a T-independent 7 ~ Cp/T when 
T > A.;; , as observed in experiments [3|] . Our estimate of 
the vison bandwidth, T^, is also consistent with a peak 
in both Cp [1] and k 01 at a temperature close to Ty. 

The vison gap, A^, obtained here is roughly the same 
as the temperature at which the l/Ti of NMR starts 
to deviate from the low temperature scaling of Eq. ^ 
0. When T is above A^,, thermally activated visons will 
proliferate. We discuss a theory of the spin dynamics in 
this thermal vison regime in the supplement, and find 
a 1/Ti with a weaker T dependence compared to that 
present for T < A„ . These observations are qualitatively 
consistent with the NMR data for 0.25 < T < 10 K 

Ref.i also measured the thermal conductivity, in an 
applied field H up to 10 T. There was little change in n for 

< 4 T. As couples to the conserved total spin, it only 
appears as an opposite "chemical potential" term for Za , 
modifying the temporal derivative {dT + {H /2)(j^)z'^ (dr-- 
{H/2)a^)z. At the quantum critical point, this term will 
induce a condensate of z i.e. a non-coUinear magnetically 
ordered state. We do not expect a significant difference 
in the thermal conductivity of the gapless spinons versus 
gapless spin-waves across this second order transition. 
We conjecture that the change at 4 T is associated with 
a vison condensation transition to a valence bond solid, 
as the field scale is or order the energy scales noted in the 
previous paragraph. This transition is possibly connected 
to the 7J-dependent broadening of the NMR spectra . 

We have described the properties of a Z2 spin liquid, 
on the verge of a transition to an magnetically ordered 
state. We have argued that the quantum critical spinons 
describe the NMR observations 12\ , while the visons (with 
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This supplement presents additional details on compu- 
tations in the main text. The large A'^ expansion of the 
nonlinear a model field theory of the transition between 
the Z2 spin liquid and the phase with non-coUinear mag- 
netic order is presented in Section I. The thermal con- 
ductivity of the spinous is considered in Section II, and 
the NMR relaxation rate at temperatures above the vison 
gap is discussed in Section III. 




I. LARGE-Af EXPANSION OF NONLINEAR a 
MODEL 



FIG. 1; Bubble diagram for the correlation function of 11. 
The solid line represents propagator given by equation (|8]). 



The phase transition between a non-coUinear Neel 
state and spin liquid state can be described by the 0(4) 
nonlinear a model as in Eq. (1) in the main text 



S=- I d^rdrUdrZj^ + C^\'^rZa\^^ 

9 



(1) 



with the constraint that |zip + |z2p = 1. The physical 
antiferromagnetic order parameter is related to the 0(4) 
field Za by the bilinear function P, 0] 

= Za^al^ai^^zp. (2) 

This differs from the collinear case, where the order pa- 
rameter is linearly proportional to the field of the 0(3) 
a model. Therefore, the spin correlation function is pro- 

Sirtional to a bubble diagram of the z field (see Fig. [T]) 

x{x,t) - T\{x,t) = (za(a;,r)z/3(a;,r)z;;(0,0)z*(0,0)) . 

(3) 

The correlation function Ii{k,ui) for the 0(4) field can 
then be calculated using the large-A^ expansion. The 
framework of the expansion can be set up in the disor- 
dered phase as follows. First, we rewrite the constraint 
as a path integral over a Lagrangian multiplier field A 

S=Ygj d^rdT[\drZa.\^+C^\VrZa.\''+i\{\Zo.?~l)]. (4) 

Here the coupling constant g is rescaled from that in 
Eq. (1) in the main text to show the N dependence in 
the large-A^ limit explicitly. Integrating out the n field 



in the above action, the path integral over A becomes 
N 



VX exp 



Trln(-c'V^-9^-f iA) 



/ drd^xX 

9 



(5) 



Therefore, in the N oo limit, the path integral is dom- 
inated by the contribution from the classical path, along 
which A becomes a constant given by the saddle point 
equation 



d^k 



1 



1 

9 



(6) 



where = iX, and m = Af"^ is the spinon gap in the 
iV — > oo limit. 

At the N — oo order, the A field is treated as a con- 
stant, and the theory contains only free Za field with 
mass gap A^. The full large- expansion is obtained 
by including fluctuations of A controlled by the action 
in Eq. ([5]): the iV~" order expansion corresponds to a 
n-loop correction. 

The spin correlation function in Eq. ([3]) can be cal- 
culated from the n(fc,w) correlation function for the Za 
field using the large-A^ expansion. At A'^ = oo order, the 
correlation function is given by a bubble diagram of two 
free propagators. 



no(A:,iw„) 



^^^^Go(p,ii^n)Go(p+fc, 



il^n+i^^n 



(7) 
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where Go{p,iuJn) is the free propagator of Za field 

At the 1/iV order, the contribution from the fluctua- 
tion of the A field needs to be included at one-loop level. 
There are two corrections that need to be included for 
the bubble diagram: the self-energy correction and the 
vertex correction. 



First, the bare propagator in Eq. ([7]) needs to be re- 
placed by a propagator with a self-energy correction at 
one-loop level [3]. 

1 

where the self-energy has two parts. The first part 
comes from an insertion of A propagator on Za propa- 
gator shown in Fig. [2l 



= ]^^Ey f^-^^ • (10) 

The second contribution is given by Fig. [31 and the total self-energy is 

11 f d/^ 

T,{k, iujn) = S(fc, iujn) ~ -g^^-Q ^ X! J ^^;;;^G'o(p,iJ^)S(fc,ii/„)Go(p,w„). (11) 

In addition to including the self-energy in the propagators of Il{k, iujn), the vertex correction (see Fig. [?]) also needs 
to be included. 

jT(iv)fi, ■ A _ 2 1 f (PpcPq Gaip, iv„)Go{p + q, iun + ien)Go{p + k, ivn + iu!n)Go[p + q + k, iv^ + + i^^n) 




FIG. 2: The first term in self-energy correction: E(fc,itj„) 
evaluated in equation (|10|l . In this diagram and Fig. |3]the 
dotted fine represents propagator of A field given in equation 
((S)), and interaction vertex between two Za field operators and 
one A field operator is given by the last term in action Q. 



A. Local susceptibility 




I 
I 
I 



FIG. 3: The second term in self-energy correction, which cor- 
responds to the last term of equation This diagram 
contains two A fields, which scales as , and also a loop of 
Za field, which contributes a factor of N . Therefore the whole 
diagram is also at the order of 1/A'^. 



In this section we consider the behavior of the imag- 
inary part of dynamical susceptibility at the threshold 
to creating two spinon excitations. Ki N — oo and zero 
temperature, the integral in the expression of Ii{k,iujn) 
can be evaluated analytically from Eq. ([7]) , and the result 



no(fc,w) 



: tan" 



_i ( y/c^k"^ — up- 
\ 2rTi 



(13) 



threshold 



1 ( uj^ ~ (?k'^ - \ir?\ , , 
Reno(A:,c^) = - In , (14) 



and 



Imno(A;, Lo) 



sgn(^) 



(cj- Vc2fc2 + 4m2). (15) 



The real and imaginary part of above equation have the 
following asymptotic behavior when lo is just above the 



Naturally, Eqs and are connected by a Kramers- 
Kronig relation. 
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FIG. 4: Vertex correction for the bubble diagram, appeared 
in equation (|12l) . 



The sharp discontinuity in the imaginary part of sus- 
ceptibiUty is an artifact of the N — oo hmit, and is mod- 
ified once we add in vertex corrections. Actually, as the 
bubble diagram evaluated in equation (|13p has a loga- 
rithmic divergence at the threshold, the ladder diagrams, 
which contains all orders of this divergence, should be 
summed as in an RPA approximation (see Fig. [5]). Since 
the bubble is divergent at the threshold, the most di- 
vergent contribution to the ladder diagrams comes from 
the propagator that is on shell and at the spinon thresh- 
old. When that happens, there is no momentum transfer 
through the A propagator. Therefore, we can approxi- 
mate the interaction vertex in the RPA summation as 
the A propagator at zero momentum: 



1 



167r 



^no(o,o) N ^^^^ 

Therefore the RPA resummation of the ladder diagrams 





FIG. 5: Ladder diagram for vertex correction. A RPA resum- 
mation of these diagrams is evaluated in equation (|16|) . 



gives 



n = Ho + nouHo 



Ho 



1 - uH, 







(17) 



Taking the imaginary part of x, and taking only the most 
divergent part, we obtain 



imn ^ 



Imnn 



u2(Reno)2' 
Pluging into Eqs ([M]) and (|15p . we obtain 



(18) 



lmU{k,uj) 



m 



(19) 



In order to relate this result for IT to the physical spin 
correlation function x, we need to insert the proportional 
constant in Eq. ([3]). Combined with the spectral weight 
from higher loop corrections, this gives the constant A 
appearing in Eq. (2) in the main text. The mass gap 
of the spinon A° = m also receives higher loop cor- 
rections and becomes in general. In addition, near 
the threshold, the factor of y/uj'^ — c^k"^ is approximately 
2/S.z- Therefore the above equation can be rearranged 
into 



Imx(A;, uj) 



16A, 



In' 



(20) 



The overall scaling Imx is a result at = oo, 

and shall be refined to Imx ~ when higher loop 

corrections are included, where fj is the scaling compo- 
nent appearing in Eq. (2) in the main text. With this 
correction, the above equation becomes 



lmx{k, Lo) 



0(\U>\ - ^C2fc2+4A2 



16A 



2-fj 



In" 



16A2 



(21) 



which is the same as Eq. (4) in the main text with C = 



B. Relaxation rate at order 1/A 

In this section we calculate the relaxation rate at 1/A^ 
order, at finite temperature above the critical point. We 
will show that, in the leading order of 1/N expansion, 
there is a singular term proportional to InN/N. 

Following Eq. (5) in the main text, we calculate 



r = lim — 



ImXL(g,w)- 



(22) 



The singularity arises from the 1 /N self-energy, by re- 
placing XL with n calculated with the full Green's func- 
tion in Eq. 

■^^-p- G(fc, Wn)G{k + q, ivn + 

(23) 

In the critical region, temperature is the only energy 
scale. Therefore we have set /3 = 1 in the above equation, 
and in the remainder of this subsection. 
Pluging Eq. ([23]) into ([22]), we obtain 



r = lim - [ '^-^^y^l-mG{k,Wn)G{k + q,Wn + i 
oj^o UJ J (27r)^ ^-^ 



(24) 
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Using the frequency summation identity 

lim — Imy^ G'i(w„)G'2(it'„ + i 
de ImGi(e)ImG2(e) 



27r 



2 sinh"' 



(25) 



and changing the variable in the second integral from q 
to k + q, we obtain 



where 



r 



A{e) 



de A{e 



27r 2 sinh 



d^k 
(2^ 



ImG(fc,e), 



(26) 



(27) 



and the Green's function includes self-energy correction 
at 1/iV order 



G{k,Lo) 



1 



I](fc,Cj) 



(28) 



Below we will see that the imaginary part of the self- 
energy leads to a InN/N term, which is more divergent 
than the 1/N correction from the real part. So if we 
ignore the real part of the self-energy for the moment, 
the imaginary part of Green's function is 



ImG(fc,tj) 



Im S(fc, Lu) 



[ImS(fc, 



r, (29) 



where E(fc,cj) is of order For the case that + 

m? — up' 7^ 0, the integrand can be expanded to the order 
of XjN by ignoring the ImS term in the denominator. 
Therefore, we expand ImS around the quasiparticle pole 
c/co — \fuP^~-Tr?: 

T.{k,U!) = S(fco,c^) 

+ J:'{ko,u;){k^-kl) + 0{{k^-k'^f): (30) 

the integral of the second term does not have a singularity 
because k'^ — k'^ is an odd function, and higher order terms 
are also not singular. Hence these terms result in regular 
corrections of the order However, integrating the 

constant term is divergent near the pole if the ImE term 
is ignored. Therefore it needs to be put back and the 
most divergent term in A{uj) is 



d^k 



ImI](/so, uj) 



(2^)2c2(fc-fco)2 + [ImE(A:o,w)]2' 
and the result of this integral is 



A(.) ^ ^ 



1 



47rc2 



arctan 



(31) 



(32) 



The function Aiuj) can be expanded to the first two or- 
ders of ImE as 



4c2 47rc2 Ci;2 _ ^2 

Pluging this into Eq. P5|) . we obtain 



(33) 



r ^ 1 



de 



1 



1 

16 



27r sinh^(e/2) 

1 ImS(\/ £2 — 777,2^ ^) 



(34) 



The first term resembles the relaxation rate in the N — 
oo limit, and the second term yields a l/(A^ln A^) correc- 
tion because the integrand diverges when e —> m 



1 



ImE(0, to) 



167r2c2 2TOsinh^(TO/2) 



In ImS(0, m). 



(35) 



Here to is the mass gap of spinon in the critical region 
a,t [3 ~ 1. In the N = oo limit it can be evaluated 
analytically 



= e = 2 In 



At order 1/iV it has been calculated that[ 



1 
r 

Thus we obtain 



lmi;(0,m) _ 0.904 



2m 



0.904 



N 



167r2c2 sinh^e/2iV 



■IniV 



(36) 



167rc2 N 
This is the \nN/N correction in Eq. (7) in the main text. 

II. SPINON THERMAL CONDUCTIVITY 

The general equation for thermal conductivity at 2-|-ld 
CFT was given in Eq. 8 in the paper: 



SC Ti, 



(37) 



The entropy density s is given in the paper by Eq. (10). 
Based on simple scaling arguments, the leading order 
scaling behavior of momentum relaxation rate l/rimp 
reads [6|] : 



(38) 



with random potential V{r) coupling to |z|2, and Vlmp 
is defined as V{r)V{r') — V^^p5'^{f— r'). For a ran- 
domly distributed impurity with Yukawa potential Vq = 



10 



Vz/{q^ + A^^) and density riimp, we can identify 



imp 



''imp 



I [I . Compensating the dimension by inserting 
the spinon bandwidth and other physical constants, 
we obtain the equation for the thermal conductivity of 
spinous (Eq. 11 in the paper): 



Nc^hk%fi''T^T, 



(39) 



Notice that this equation is only applicable to the case 

with 1/Timp <C T. 

III. THERMAL PROLIFERATION OF VISONS 

In this section we discuss the regime T > A^, where 
the visons have thermally proliferated. As noted in the 
paper, at these temperatures the 1/Ti NMR relaxation 
rate is observed to have a plateau j5| . We believe this is 
a general feature of a dense vison regime: the presence 
of visons makes it harder for the spinous to propagate 
independently, and so a vector spin model (which has a 
T- independent NMR relaxation rate j^]) becomes more 
appropriate. 

Here we will illustrate this qualitative idea in a spe- 
cific model. Rather than thinking about this as high T 
regime for visons, imagine we reach this regime by send- 
ing A„ ^ at fixed T. In other words, we are in the 
quantum critical region of a critical point where the vi- 
son gap vanishes leading to phase with the visons con- 
densed. We have already argued in the paper that the 
spinous are also in the quantum critical region of a tran- 
sition where the spinous condense. Thus a description 
of the spin dynamics in the regime T > A^, is provided 
by the quantum criticality of a multicritical point where 
both the spinous and visons condense. A general theory 
of such multicritical points has been discussed in a recent 
work by two of us. 2] The NMR relaxation is then given 
by 1/Ti ^ yjmc^ where r]mc is the anomalous dimen- 
sion of the magnetic order parameter at the spinon- vison 



multicritical point. 

Our only present estimates of r]mc come from the 
1/iV expansion, and so it is useful to compare estimates 
of anomalous dimensions in this expansion at different 
quantum critical points. For the regime, T < A^, dis- 
cussed in the main paper, the NMR relaxation is con- 
trolled by the theory in Eq. ([T]) describing the condensa- 
tion of the spinous alone. Here we have l/Ti ~ where 
i 



1 



64 



(40) 



In the higher temperature regime, T > A^, we have 
l/Ti ^ T^'fmc^ and the same 1/N expansion for this ex- 
ponent at the multicritical point where both spinous and 
visons condense yields 



256 



iir^N l + 256/cV(7r2Ar2)- 



(41) 



Here k is the level of the Chern-Simons theory describ- 
ing the multicritical point. The large N expansion is 
performed with k proportaional to N, and the physical 
values are k — 2 and iV = 4. 

The key point is that rjmc < rj. Hence 1/Ti will have 
a weaker dependence on T for T > Az than for T < A^ . 
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